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CONTROL OF SOLID-APOGEE MOTOR DURING ORBIT TRANSFER MANEUVER 


Zhang Yuntong 

(Beijing Institute of Spacecraft System Engineering) 


Abstract 

This paper provides a method for determining the optkr.’.ui time and 
attitude in process of geostationary satellite launch at which the 
solid-apoaee motor to be fired. 

Subject Terms: firing time, satellite attitude, optimum design. 

Apogee engine 

(The manuscriDt was received on July 13. 1989.) 

I. Forward 

In the process of launching geostationary satellites, the Apogee 
engine on the satellite provides an increment in the velocity of the 
satellite which, ultimately, transforms the orbit the satellite 
travels from a transitional orbit into a geostationary orbit. The 
direction of spin of a geostationary satellite is the direction of the 
velocity increment. Therefore, the orbit transition control of an 
Apogee engine is to determine the firing time and the Spin direction 
(the so called firing attitude) for a given transitional orbit. 

Different firing times and firing attitudes can result in different 
orbits after the engine is stopped. Hence, the firing time and firing 
attitude should be carefully chosen to obtain the optimum orbit. In 
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the actual engineering practices, the choice is always subjected to 
some limitations, and. therefore, this problem can be thought of as an 
constrained optimization problem. In this paper, we will introduce: 

(1) the working model of Apogee engines, (2) the definition of 
optimization target functions, and (3) method of optimization. 

II. The Working Model of Apogee Engines 

1. Equation of Motion 

The equation of motion of the center-of-mass of a satellite during 
the working period of Apogee engine can be simplified as: 


d 1 * - _ + F o 

dt* r* + “ 

d*y _ n ■ F Agtotil 


where x. y. z are the coordinates of the satellite in the earth 
equatorial coordinate system: J* is the distance from the satellite to 
earth's center: * is earth's gravitational constant: F is the thrust 
of the Apogee engine; m^ and m are the initial mass of satellite 
and the rate of mass flow of propellent (per second); j tiihd are 
the longitude and altitude of the firing attftutts. 3he 
should be carried out in the time interval of the working period of 
the Apogee engine; namely, the integration starts at the ♦’-♦SljfvWhen the 
Apogee engine is ignited and stops when the Apogee engine is stopped. 
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After the transformed orbit is selected for a firing time 


the 


initial conditions for the set of equations (1) can be obtained and F. 
mg. m and the working time of the Apogee engine are all fixed 
parameters. If a set of.• and i were chosen, the orbit of the 
satellite after the engine is stopped can be obtained from the 
integration of equation (1). In other words, for a set of firing time 
and firing attitude (tg ,«»£'), the orbit of the satellite after the 
engine is stopped can be calculated. 

2. Pulsed Orbit-Transition Model 

For solid-Apogee motors, the working time is on the order of a few 
tens of seconds. During*its working period, the satellite is in the 
vicinity of the Apogee point. It can be assumed that the satellite 
obtains the velocity increment instantly. The position and velocity of 
the satellite at ignition and stop of the Apogee motor should be 
expressed in terms of the subscripts "0" and "f". Their relationship 
can be simplified as 

*. = */ 

»€=*/ 

*,=^*+FoaB3*inc 

where V is the velocity increment provided by the Apogee motor and it 
quantity is.- 

• ■ ■ 

where g 0 is the gravitational acceleration on the surface of the 
e a rth: I sp is the specific impact in vacuum; and T is the working time 
of the motor. Similarly, if the firing time tg is determined. x 0 . yg, 
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zq, xq. yO' ZQ can be calculated from orbit measurement and a and 6 
are the quantities which should be chosen. The relationship between 
the parameters x s . y s , z s . x s , y s . z s and the orbital parameters is 
determined by the well known equations as those in [1]. 

3. Comparison and Conclusion of The Two Models 

In the two methods mentioned above, method (1) is closer to the 
actual motion. When the working time of solid motors is a few tens of 
seconds, the difference of the calculated target functions can be 
ignored. The minimum difference in the target function is obtained 
when the difference in the firing times based on these two models is 
0.55 times of the working time of the motor. Therefore, optimal 
condition can be obtained if the actual firing time is 0.55T in 
advance of the time based on the pulsed model. 

III. Tai'cet Function 

1. Definition of the Target Function 

In order to evaluate the advantages and disadvantages of the chosen 
parameters tg. a • and 6, the optimal target function is defined as: 
the corrective velocity increment necessary to bring the orbit at 
engine stop to the geostationary ncbit, The optimum setof functional 
variables is the set of tg. a, and 5 Stfueh minimizes the target 
function. ' 

2. Calculation of the Turgct Function 
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The orbital parameters at engine-stop are: geometrical 1onaitode a 
enaine-ctoo A,, halt 1 ona-axie of the orbit «,i orbital eccentric 
ratio e g : oibital tilt anqle j g : equatorial longitude of the up-r.tsin 
point Q.. The parameters of the geostationary orbit are: fixed 
position Geometrical lonaitude ' half long-axis of the orbit * 
orbital tilt angle i Q : equatorial longitude of the up-rising poim 
In order to maintain the north-south positionshould be around 
270° and i 0 should be about 0.1°. To maintain the east-west position. 
e s should be less than a certain value and this value was taken a? 
zero when the target function was calculated. The corrective velocity 
in the target, function was expressed in terms of three parts: 

(1) Corrective velocity AK, necessary for the correction of half 
lone-axis and eccentric ratio: 

To correct a,, e s to e..(eo“0), a dual-pulse was required and the 
veiocitv was increased (or decreased) at the hunch point. The 
caIcu1 atlona1 equation? are: 


r. p = «.(! — «,) 



The corrective velocity + W ,—FJ 

(2) Corrective velocity A V 2 iwcessary for fixed-point capture: 
Assuming that and in order to maintain the satellite at 

position , the satellite should drift towards the west after the 
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motor was stopped. If the rated drifting velocity is 4(degrees/ day ), 
1<0 . The half long-axis corresponding to the drift should be 




The ca1cu1 ationa1 equations for corrective velocity AFj according to 
various a. and e s are tabulated in table 1: 


Table 1: Calculational equation for corrective velocities 

mrnmm**** 



AT, 

«.>a< 

0 

«.<<>< r„>a, j 

•>Z (04-«.) 

r, p <a 0 I «,>04 

0 

! | 

•l «««'•,> 

! «i<*. I 

•z (Mi - «,) 


'kev 1: condition’ 

y 

Ir: table 1. w« = { is the anaular spet-d of earth's SDin and 

(3) Corrective velocity Af, necessary for orbital plane capture: 
Orbital plane capture is achieved by correcting ■§ to io and4?. i 

and the orbital plane should be tilted an i D angle. 

The corrective velocity can be appr«*MM* *«: 
when i p < ir/2. AV 3 -V()ip 

when ip >*/2. ^ V 3 “Vq (ir—i p ) (i p in radians) 

If A s f . Aq • similar results can also be obtained. 

* 




3. Target Functional Values 

The target function is A V- AVj + AV 2 * *^3 an d the relationship 
between the chosen parameters (to,“ . 4 ) and the target function is 
explained in the above sections. 

IV. Optimization Process 

1. Method of Optimization 

Under the normal conditions, in the process of satellite iaunchina 
the error encountered m the transition orbit is usually verv small 
and the reaion for opt miizat ion is also a very narrow reaion: on th*- 
other hand, the calculation based on the pulsed model is also very 
simple. Based on these reasons, the method of optimization can be the 
reliable listincr method and the searching can be carried out in the 
vicin 1 tv of the standard values. 

2. Calculation of the Standard Values 

It was assumed that the standard values satisfy the two conditions 
1 ) 1 s . fi g equal i Q , respectively; 2 ) end V g are mutually 

perpendicular and x s x g +y s y s +z g z s -0. Further, it was assumed that the 
instantaneous roots at the uD-rising point before ignition. the 
motor were 03 . 03 , 13 , w 3 , # 3 , H 3 , f*#placing the instantaneous roots 
of at the ignition point. The caiculefcSen of the etendard value* of 
the firing time and firing attitude was then based on the familiar 
two-body equations: 
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oos», = oo«r 0 oo8»»+ ani,mni.amW,~ Q t ) 
*n«„= an(£,-fl,)dni,/dni. 


*= 2 ' ig "[v&(l)] (I»I "•*)■© 

Af = Af, + E,-e,sinE, 


p,=«,d-«i) 

r,«P,/(l + f i C08/,) 



e,= «m-‘Cv/ MPiA.tY! (fi>* ■* *. *<3D 

a,= sm _1 ( sin», sin («„+«*)) 

•i=Ci+ tg'‘(00Bi»tg(a,,+e,)) Ctg' 1 WUH-^ + 
008?,= sine, cos*, 

008 X= -Ctg£,ctg?, 

sin?, = p tin?, 
cot n= -tga,clg(«j«+£») 

an 6 = «infl,cos(?, + y,) + oo8A,rti(?, + ?,)cat(r?+x) 
sin (a,-7) = sin (?, + ?,)(*! Or+*)/oo8« 


key 


1: £ 3/2 and E 3/2 are in the same quadrant. 

2: When f3 is greater than » Ue 3 lies in the second quadrant. 
3: tg -i and U3 q+ £ 3 lie in the ttmw quadrant. 


3. Conclusion of the Optimization Process 

The control of the orbit transfer maneuver of Apogee motor can be 

concluded as: 

( 1 ) measurement of the transition orbit: 







< 2 ) calculation of the standard values r, fc, 7, 

(?' confirmation of the region-of-choice and step-length of t.a . 
and 6 : 

(4) calculation of the parameters a s - « 8 . i # . after stop of th 
motor for each set of t, q^and 6 ; 

(5) calculation of the target function A V; 

( 6 ) repeat of steps (3) ^(5) to determine the t, o , and 6 
corresponding to the minimum A V; 

(7) calculation of t-0.55T and (a.. & ) for actual firing time and 
firing attitude, respectively. 


REFERENCE? 
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APPLICATION OF SPLIT-PHASE COPE P3K SIGNAL IN DATA TRANSMISSION OF 

METEOROLOGICAL SATELLITE REMOTE SENSING IMAGE 

Dai Ying 

(Shanghai Institute of Spacecraft Engineering) 

Abstract 

In this paper, split-phase code PSK signal is analyzed, its 
character is reviewed, the basic formation method is briefly 
described. The application of the signal in digital image transmission 
of astronautical remote sensing is introduced. And finally, photos of 
real-time data transmission are presented. 

Subject Terms: phase shift keying, phase modulator, picture 

transmission, meteorological satellite, application 

(The manuscript was received on May 5. 1989.) 

(The definition of split-phase code can be found in reference 1.) 

I. Split-Phase Code Phase-Shift Keying Signal 

Split-phase code phase-shift keying (PSK) signal is the modulated 
signal of the carrier wave, based on the split-phase code signal, 
after phase-shift keying. It is a binary modulated signal; namely, the 
PCM-PSK/digital split-phase, which is also a primary phase-modulated 
signal without the auxiliary carrier waves. Its mathematical formula 
is : 

SO) «=•.«■ &*+***«>} (1) 

K> 






where u m - amplitude.- 

<*>*- - frequency of the carrier wave ; 

- split-phase code phase-shift amplitude. 
c(t) - split-phase code, its value should be taken as +1 or -1 
randomly. 

The relationship between the amplitude and phase-shift angle of the 
split-phase code PSK signal is: 

JJ*? ( 2 ) 

where u cm is the voltage amplitude of the modulated carrier wave. 

In the process of product research and development, according to 
equation (2). the phase-shift amplitude of the carrier wave by the 
random solit-phase code can be measured with the frequency 
spectroscope. Using this method, the modulation of the transmitter can 
be measured with higher accuracy. This equation can also be used in 
the orbital test of satellites. 

The mathematical expression of the frequency spectrum of split- 
phase code PSK signal is: 

F(a>-*> c ) = *a» l A*^$(w-«s) 



where I is a function with screening characteristics which is also 
called the 5-function and T is the period of the split-phase code. 

According to equation (3), the frequency spectrum can be plotted as 
shown in figure 1. The shape of coverage of the spectrum of this code 
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is similar to that of the fundamental band signal frequency spectrum 
of the solit-phase code. Therefore, in order to avoid interference, it 
was only necessarv to filter the visual wave, before modulation, to 
accomplish the intra-band constraint of the spectrum of split-phase 
code PSK signal. Technically, this method is easier to perform than 
the intra-band constraint of the microwave after modulation and. also, 
this method is more economically feasible. 

The split-phase code PSK signal has certain characteristics which 
makes it very powerful in the digital image transmission. Firstly, 
since the phase shift angle of the split-phase code was not equal to 
"n/2. the carrier wave part could be saved after modulation. Secondly, 
since the part of the fundamental frequency spectrum of the split- 
phase code at zero frequency was zero and the phase-shift modulation 
of the carrier wave was the linear transformation, the part of the 
carrier wave after modulation existed independently and the 
interference of the chanae of the messaae spectrum did not exist. This 
led to the features of stable carrier frequency part and the high 
signal to noise ratio. For example: when the narrow band in the 
vicinity of the carrier frequency in the frequency spectrum of the 
split-phase code A/=o.35^-tt, the power of the message part of the 
continuous spectrum was less than 0.03% of the total power of the 
message spectrum and the actual signal to noise ratio could be as high 
as 75dB. 

In figure 1. the percentages were the results of the numerical 
integration of equation (3) and indicated the power of each specific 
wave peak as compared to the total power of all the wave peaks in the 
continuous wave spectrum. One can see that most of the power of the 
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ontinuous frequency spectrum 
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1: Frequency spectrum of the split-phase code PSK 






continuous wave spectrum existed in the first wave peak which amounted 
to 90%. 

On the other hand, under the condition of total symmetry of the 
split-phase code, there was no dispersion frequencies of the modulated 
frequency spectrum of the split-phase code PSK signal. Also, there was 

no even parts. In this way. the erroneous tracking of the carrier wave 

. 

by the ground station could be avoided. 

With these special features, the split-phase code P9K signal was a 
useful tool in determining the modulation functions of a transmitter 
and the design of the match of system frequency band. If the band 
width which occupied 90% of the energy was 2.8MHz, considering 1) the 
frequency deviations caused by the carrier wave and the signal itself, 
2) the acceptable signal to noise ratio and, 3) the deformation of the 
waveform of the signal, the band width of the receiver was chosen to 
be 3MHz. Even though the edge of this band width overlapped with the 
corresponding frequency spectrum of the TIROS-N satellite of USA, no 
interference should occur because of the separation in space and time. 

If the special features of signal mentioned above were linked with 
the system function, the advantage of the split-phase code PSK signal 
was made even more obvious. For example, to transmit digital image 
with split-phase code PSK signal, when the receiving/demodulat ion 
equipment received the PSK signal, the PSK signal power which entered 
into the carrier wave lock-up phase loop of the receiver was very 
limited. In other words, even though the signal received was the 
modulated signal, the signal to noise ratio of the part of the carrier 
wave which entered into the loop was high. In this way, the signal 
could be transmitted under the condition without the modulation of 






auxiliarv carrier wave and the receiver could perform tracking of the 
carrier wave and de-modulation of the signal simultaneously. 

The disadvantages of the split-phase code PSK signal was that when 
the signal was transmitted through wireless transmission, the 
interference would occur due to the effect of multi-path electrical 
transmission. Therefore, when the system was tested or measured, the 
geographic effect should be adequately taken into account. 

II. The Formation of the Split-Phase Code PSK Signal 

The split-phase code PSK signal was rormed oy Keying the phase of 
the carrier wave signal with the split-phase code. In the following, 
two basic methods for the formation of split-phase code PSK signal 
would be bneflv explained and mathematical methods were used to 
analyze the formation process. 

Method 1: 

The basic principle is shown in figure 2 and the relationship 
between various phases is shown in figure 3. 

Figure 2 shows that the multiplier does not shift the phase of the 
carrier wave signal u(t) in the first half of the cycle of the split- 
phase code operation time, but in the second half of the operating 

cycle, u(t) is shifted 180 degrees. The mathematical analysis is: Let 
T 

u (t) “UfneosOot, when the output signal of the -multiplier is 

•»<*)=t*.cce<<M 

where k is a constant. 

At the same time, the output signal of the adder is 
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shifter 


split-phase code signal 
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sol it-phase code phase shifting principle 


HjH) " 

n< 


Fiaure 3: Phase relationship of various components 





S(t)xN.oap«M-«>,rin&>|! 

* .4.008 A<p) 


(4) 


where 


A m =%/*£♦ **• t 

4 * =•"** tefl* 


Therefore, the phase-shifting modulation of the carrier wave by the 
positive electrical potential (subei&iWy code “l") of the split-phase 
code wos accomplished. When y<*<7V , u(t) was shifted 180 degrees. 

y +f) 

j(0=«(o+a,<0 

= ».0aM%f+ir t »4natf 
= .4 .««<«,*-A<p) 


where » 

- Af>=arctg J^|. 

Therefore, the phase-shifting modulation of the carrier wave by the 
negative electrical potential (subsidiary code "0") of the split-phase 
code was accomplished. It can also be seen that by controlling the 
accuracy in phase-shifting and the amplitude ratio of the two signals, 
the split-phase code PSK signal can be formed. 

Method 2: 

Figure 4 shows a reflective-type modulator which aceooplishes the 
phase keying control of the carrier wave by controlling tbe path L of 
the carrier wave. When the split-phase code was at the positive 
electrical potential. the PIN tube was connected and point A was 
arounded and the carrier path was Lj. When the split-phase code was at 
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micro band matching network 
split-phase code 



Fiaure 4: A reflective phase modulator 



the negative electrical potential, the PIN tube was disconnected and 
the path of the carrier wave was L^- If •^5 (Lj-ly-Af is sat isfiei , 
the carrier phase transformed from f>,to fitf can be obtained. IP* was the 
initial value. 

In actual application, the 3dB bridge was used in replace of the 
cyclotron in figure 4. In fact, the modulator manufactured in this way 
would perform better and actual manufacturing was accomplished. 

III. The Application of the Split-Phase Code PSK Signal 

The split-phase code PSK signal was used not only in the digital 
communication of spacecraft, it was also used in the transmission of 
high resolution remote-sensing images of TIROS-N meteorological 
satellite. The meteorolocrical satellites manufactured in China also 
used the split-ohase code PSK sianal to transmit ultra-high resolution 
cloud maps. This system consisted of two parts (satellite transmission 
a nd or :■ und r e :• eiving » : 

Satellite Transmission: The satellite transmission included the 
electronic circuitries of the remote sensing system and the 
transmitter and antenna of the transmitting system. The remote-sensing 
system installed on the satellite was a five-channel scanning 
radiation monitor which received the thermal radiation data from the 
emission and reflection of the earth-atmosphere system and the visible 
light data which included the meteorological data from cloud top, 
earth surface and ocean. 

The remote-sensing equipment first transformed the gathered data 
into electronic signals and then quantized and digitally processed 
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these sighais. After coupling the service code, the information was 
edited and according to the special format requirement, such as the 
service code of synchronous signal, satellite attitude, earth 
meteorological data, and the error message detecting sequence, the 
binary code was transformed into the split-phase code. After phase- 
shifting keying of the carrier frequency of the transmitted cloud map 
data, the split-phase code PSK signal was generated. This signal was 
transmitted to ground, via antenna, after amplification. 

Ground Receiving: After the ground station received the split-phase 
code PSK signal, the high-frequency part was amplified, the 
frequencies were mixed, intermediately amplified, and then transmitted 
to the lock-phase de-modu 1ator. During de-modulation, the phase of the 
£■; irial was compared with the phase of the split-phase code PSK signal 
carrier wave reference signal and the split-phase code signal was 
obtained . The ,-.>::p 1 1 f led split-phase code was transmitted to the bit 
svnchror. iner to extract the synchronous clock signal and the spl i t - 
phase ie was transformed into the non-zero code and transmitted to 

the frame synchronizer. In the frame synchronizer, the beginning 
of a frame sianal—the 60 bit pseudo-random signal was examined to 
extract the frame synchronous timing signal and the cloud maps or the 
earth meteorological signal was separated into five channels. For the 
error message signal at the end of a frame of signals, the error rate 
would be examined and rate of error message of the whole string of 
data would be derived. If the error rate satisfied the requirement, 
the cloud map signal would be transmitted to the data processing 
system for real-time image processing. Or, the data could be stored 
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for future processing to obtain the meteoro)oaicel data of the ci-un 

top. earth surface and ocean. 

Figure? b and 6 show two processed high resolution cloud maps which 

were transmitted by the high resolution cloud map transmission system 

1 HRPT svsterro of the "Fong-Yun" No. 3 meteorological satellite via 

split-phase code PSK signal. Figure 5 shows the cloud map transmitted 

bv the HRPT system at 08GMT on September 8. 1988. It shows the rain- 

producing cloud system of the south eastern part of China. In the 

southern part of China, there were several cloud strings which fed to 

the rain-producing cloud region. The dark region at the bottom of the 

figure was the South Sea and Hainan Island was at the center of it. 

Figure 6 shew? the visible light high resolution meteorological cloud 

man transmitted bv HRPT svstem. at the downstream of the Yangtze 

river at 6:30 GMT on September 20. 1988. This map shows clearly the 

o! ovd -listrib-it ion and the sand content of the downstream of the 

Yangtze rive and the East Sea. It also shows the position of the five 

' o’; .’ce &*>■•! the width of each string at the coastal region of the 

diangSu province. Also, the south pointed dispersion of the fresh 

water from the Yangtze river in the ocean can be clearlv seen. This 

dispersion of the fresh water affected the HangJu bay. Further, the 

mixing and distribution of the fresh water and sea water can be seen. 

(Figures 5 and 6 are in the third cover page.) 

All of these figures have clear images and clear geographic 

boundaries. The resolution was less than 1km. One of the high 

r 

resolution cloud maos was used by the "AVIApN WEEK” magazine. 

Data transmission was one of the key issues to produce high 
quality cloud maps. Therefore, it should be noted that transmission of 
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the remote-sensing imaces by the split-phase code PSK signal was 
successful. These cloud maps provide important data for weather 
forecasting, typhoon monitoring, ocean detecting, and agricultural 
evaluation. 
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